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PREFACE 


This  is  the  second  of  a  two  volume  report  on  the  subject  of  friction  damping 
of  dynamic  systems.  The  first  of  these  reports  entitled  "Study  of  Character¬ 
istics  of  Dry  Friction  Damping"  was  prepared  as  an  annual  report  in  March  1984 
and  submitted  to  the  Air  Force  Office  of  Scientific  Research. 


The  first  report  contained  an  introduction  to  the  topic  along  with  a 
thorough  review  of  pertinent  literature  and  discussion  of  (a)  interface  charac¬ 
teristics,  (b)  characterization  of  friction  forces,  (c)  coefficient  of  friction 
and  (d)  dry  friction  damping  technology.  An  analytical  model  and  a  corresponding 
experimental  set  up  was  also  presented  as  the  basis  for  further  research. 

Possible  modifications  to  a  nonlocal  lav  were  also  included.  It  is  recommended 
that  the  reader  study  the  first  report  prior  to  commencing  the  present  volume, 
although  this  report  is  essentially  self  contained. 
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SUMMARY 


This  report  contains  results  of  analytical  and  experimental  efforts  related 
to  dry  friction  damping.  A  general  time  history  solution  is  developed  for  a 
multi-degree-of- freedom  system  which  incorporates  the  use  of  nonlocal  lavs  of 
friction.  The  formulation  and  the  solution  are  entirely  nonlinear  and  have  the 
necessary  flexibility  for  modelling  specific  engineering  problems  where  the 
presence  of  friction  at  interfaces  influences  the  dynamics.  Two  examples  are 
chosen  to  illustrate  certain  important  aspects  of  the  procedures. 

A  friction  test  assembly  designed  and  fabricated  under  this  program  served 
as  an  experimental  rig.  The  results  from  a  limited  number  of  tests  are  presented 
and  discussed.  A  comparison  between  analytical  and  experimental  results  is 
given. 


This  annual  and  final  report  covers  the  period  between  March  1,  1984  and 
April  30,  1985.  The  work  reported  here  was  under  the  sponsorship  of  the  AFOSR 
under  the  technical  direction  of  Dr.  A.  Amos. 
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1 .  INTRODUCTION 


One  of  the  most  important  sources  of  energy  dissipation  in  built-up  struc¬ 
tures  pertains  to  dry  friction  at  joints  and  at  interfaces  in  contact  with  each 
other.  When  materials  with  very  low  inherent  damping  are  used  in  environments 
offering  essentially  no  aerodynamic  damping  (such  as,  for  example,  in  space 
structural  systems),  the  ability  of  the  structures  to  withstand  resonant  vibra¬ 
tion  depends  almost  entirely  on  the  extent  of  friction  damping.  Even  in  the 
design  of  components  which  can  generate  certain  levels  of  aerodynamic  damping,  it 
is  essential  to  be  able  to  estimate  accurately  the  levels  of  non-aerodynamic 
sources  of  damping  because  unlike  the  latter,  the  former  can  be  either  positive 
or  negative.  For  example,  the  aeroelastic  instabilities  of  aircraft  structures 
such  as  airplane  wings  or  turbomachinery  blades  are  attributed  to  negative  aero¬ 
dynamic  damping  developing  in  the  system  as  a  result  of  its  interaction  with  the 
air  forces  causing  the  system  to  vibrate.  In  such  situations,  escalation  of 
vibratory  stresses  to  dangerous  levels  is  prevented  only  through  the  contribution 
to  damping  from  sources  such  as  friction.  In  some  instances,  artificial  devices 
are  designed  and  introduced  into  the  structural  system  with  the  sole  purpose  of 
developing  friction  forces  leading  to  damping  when  the  dynamics  of  the  component 
demands  it.  Thus,  there  is  a  clear  need  to  enhance  the  level  of  understanding  in 
the  broad  area  of  dynamics  of  surfaces  in  contact. 

Among  the  multitude  of  factors  contributing  to  the  dynamics  of  elastic 
bodies  in  relative  motion  that  come  into  contact  with  each  other  at  interfaces, 
the  most  significant  factor  is  the  governing  law  of  friction.  It  is  the  law  of 
friction  that  accounts  for  the  characteristics  of  contacting  surfaces  and  the 
environment;  and  establishes  the  relationship  among  forces  and  resulting  motion. 
The  success  of  prediction  of  dynamic  response  of  a  structural  component  depends, 
to  a  large  extent,  on  the  choice  of  a  realistic  friction  law  that  accurately 
portrays  the  interfaces.  The  laws  of  friction  so  constructed  must  be  adaptable 
in  a  dynamic  analysis  of  the  structural  system  and  directly  verifiable  in  a  care¬ 
fully  controlled  test. 

This  report  is  aimed  at  examining  different  characterizations  of  friction 
forces  at  interfaces.  A  general  nonlinear  theory  of  dynamic  systems  is  developed 
in  such  a  way  that  a  local  or  nonlocal  law  guided  by  experimental  data  may  be 
incorporated.  Solution  procedures  are  presented  and  discussed. 

This  necessarily  includes  an  extension  of  simple  "local"  characterization  of 
friction  for  the  neighborhood  of  a  single  surface  asperity  to  a  "nonlocal"  or 
gross  characterization  for  an  interface  having  a  large  number  of  asperities. 

Three  laws  of  friction  that  are  adaptable  to  dynamic  analysis  of  structural 
systems  are  outlined  and  discussed  below. 
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2.  LAWS  OF  FRICTION 


One  of  the  laws  considered  here  is  the  simplest  and  represents  the  most 
commonly  used  form  i.e., 

^  • 

F  -  F  sgn(X) 


This  law  is  used  when  |xj  is  above  a  specified  threshold  of  "sticking"  and 
is  replaced  by  a  "static"  friction  condition  if  the  relative  velocity  is  below 
the  threshold  (this  will  be  discussed  further  below). 


A.  Friction  Force  of  Constant  Amplitude 

An  application  of  this  law  for  the  system  shown  in  Fig.  1  served  as  the 
basis  for  the  parametric  study  leading  to  the  design  of  an  experimental  rig.  Tt 
latter  will  be  discussed  in  a  subsequent  section  but  the  highlights  of  the 
analysis  are  presented  below. 

The  governing  equation  of  motion  for  the  system  shown  in  Fig.  1  may  be 
written  as 


Mx  +  K(x-X)  +  F  sgn(x)  +  C(x-X)  «  0 


(2.1 


which  may  be  written  in  a  nondimens ional  form  as 


2«q  ♦  CZq 


-  sgn(q)D  -  -W  f (t) 
C 


(2.21 


where 


Mid2  _  _ 

q  *  -  x,  x  **  x-X,  X  ■  XQf ( cot ) ,  T  *  wt 

**0 


The  derivatives  in  Eq.  2.2  are  with  respect  to  t. 


,«*  ft' 


R85-956479-2 


The  numerical  solution  of  the  governing  differential  equation  2  is  obtained 
using  the  Newmark  integration  scheme  (Ref.  1)  with  discrete  time  steps.  The 
calculations  are  performed  when  the  mass  M  is  moving  under  the  action  of  either  a 
sinusoidal  or  step  input  displacement  function.  When  the  block  is  stuck  (as 
indicated  by  the  relative  velocity  dropping  below  a  pre-specif ied  threshold) 
static  friction  conditions  are  imposed.  A  test  is  then  made  at  each  time  step  to 
determine  if  the  accelerating  force  on  the  block  exceeds  the  static  frictional 
force.  As  soon  as  this  occurs,  dynamic  conditions  are  imposed  and  the  motion 
parameters  are  calculated.  A  check  on  the  magnitude  of  the  velocity  indicates 
whether  or  not  the  block  has  become  unstuck.  If  the  velocity  threshold  has  not 
been  exceeded,  then  static  friction  conditions  are  reimposed.  This  process  is 
repeated  until  the  block  slips. 

Two  points  are  worthy  of  note  in  this  type  of  solution.  The  first  point 
pertains  to  the  definition  of  the  "stuck"  condition.  With  the  present  approach, 
the  friction  force  alternates  as  the  velocity  bounces  within  the  specified 
threshold  (±)  region  causing  spurious  acceleration  spikes.  The  size  of  the 
velocity  threshold  has  a  marked  effect  on  the  response  and  guidelines  may  need  to 
be  obtained  from  test  data.  The  second  point  is  that  the  model  permits  incor¬ 
poration  of  different  friction  models  such  as,  for  example,  microslip  in  place  of 
the  "stuck"  condition. 

As  stated  earlier,  the  model  was  used  in  the  design  of  the  experimental 
setup.  Figures  2  and  3  show  typical  results  obtained. 

These  results  were  obtained  using  the  parameters  for  the  experimental  rig 
i.e.  the  mass,  spring,  normal  load,  applied  g  loads  etc.  correspond  to  those 
used  in  the  test.  Figure  2  shows  the  input  displacement  and  the  resulting  dis¬ 
placement,  velocity  and  phase  plane  plot.  The  flat  top  in  the  output  displace¬ 
ment  characteristics,  the  nature  of  stick-slip  motion  evident  in  the  velocity 
characteristic  all  lead  to  the  stable  limit  cycle  shown  in  Fig.  2.  The  kinks  in 
the  phase  plane  plot  represent  the  displacements  and  velocities  at  which  the 
system  sticks  and  slips.  Figure  3  shows  the  variation  of  friction  force  as  a 
function  of  time,  relative  and  absolute  accelerations.  The  latter  will  be  used 
in  a  later  section  of  this  report  for  making  comparisons  with  measurement. 


■r.ss.y.-r 


^'.NV  . 
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B.  Friction  Force  Increasing  with  Relative  Velocity 


Another  law  that  was  examined  in  this  study  is  of  the  form 


V, 

F  ■  F^tanh  — *•  (2.4) 

° 


This  law  has  the  advantage  in  that  it  degenerates  into  the  simple  law  presented 
above  or  the  familiar  viscous  law  with  appropriate  assumptions.  The  single- 
degree-of- freedom  system  shown  in  Fig.  4,  differs  from  that  of  Fig.  1  in  that  the 
spring  and  viscous  damping  forces  relate  to  absolute  motion  of  the  masses  rather 
than  to  motion  relative  to  the  frame.  The  equation  of  motion 


Mx  +  kx  +  cx  *  F  with  Vr  M  a  -  x, 


(2.5) 


can  be  solved  numerically  using  a  simple  Euler  integration  scheme.  Assuming 

Fq  =  0.1  Ka0 

0)  *  0.1  (D0,  u>0  *  /k/M 

VQ  “  0.01  a0% 

C/CQ  *  0.02 


the  computed  response  is  shown  in  Fig.  5  for  a  typical  cycle,  indicating  that 
the  assumed  friction  law  results  in  a  high  frequency  ringing.  An  attempt  was 
made  to  model  an  experimental  set-up  shown  in  Figs.  6  and  7,  because  test  data 
from  this  set-up  had  shown  responses  with  high  frequency  ringing  (see  Fig.  8). 
This  had  been  attributed  earlier  to  possible  electronic  noise.  However,  with  the 
hyberbolic  tangent  law  it  was  possible  to  reproduce  these  responses  by  modeling 
the  experimental  hardware  as  shown  in  Fig.  9.  The  elasticity  in  the  dynamic 
system  is  allowed  in  the  model.  The  equations  governing  the  motion  of  this 
system  are 


•  •  * 

mXj  +  cXj  +  kX1  ■  F 
•  •  * 

mX,  +  cXo  +  kX9  =  P.sin  wt  -  F 


(2.6) 
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Again,  the  friction  law  was  chosen  to  be 


•  • 

X  -X, 

F  *  F„  tanh  -2 — i 
°  Vo 


(2.7) 


These  equations  were  integrated  using  the  Euler  integration  scheme  and  typical 
time  histories  of  response  are  shown  in  Figs.  10,  11,  and  12  for  various  forcing 
frequencies.  The  viscous  damping,  was  chosen,  as  before,  to  be  2  percent  criti¬ 
cal  with  FQ  =  0.1  PQ,  VQ  =  0.01  PQ//Km  and  Q  *  u>/ii>0.  The  ringing  is  clearly 
evident  in  the  response  coordinate  .  The  X2  coordinate  is  a  smooth  sinusoid  as 
it  is  being  driven  by  the  forcing  function.  As  the  forcing  frequency  approaches 
resonance,  the  ringing  disappears  but  relative  motion  at  the  interface  persists, 
resulting  in  damping.  The  results  from  this  study  show  that  modeling  which  takes 
into  account  the  elasticity  of  the  system  may  produce  responses  that  are  more 
realistic.  Thus,  in  reducing  data  from  complex  structures  care  must  be  exercised 
so  that  certain  responses  may  be  recognized  as  inherent  to  the  structural  system 
rather  than  mere  noise  that  should  be  filtered  away. 


C.  Extension  of  Nonlocal  Friction  Law 

Finally,  the  nonlocal  law  postulated  by  Oden  and  Pires  (Ref.  2)  and 
discussed  in  the  first  volume  of  this  report  (Ref.  3)  has  been  extended  to 
include  velocity  dependence  for  the  forces  at  the  interface  and  to  integrate 
along  the  path  followed  by  the  points  at  the  interfaces  undergoing  relative 
motion.  Thus,  the  velocity  dependence  accounts  for  the  dynamics  of  relative 
motion  at  interfaces  and  the  type  of  integration  proposed  recognizes  a  noniso¬ 
tropic  weighting  function.  The  nature  of  this  extension  is  presented  below. 

Consider  two  surfaces,  A  and  B,  in  relative  motion  and  assume  that  surface  B 
is  stationary.  The  two  surfaces  may  have  many  points  in  contact  with  each  other 
at  time  tj .  Let  Xj  be  the  surface  coordinates  on  B.  The  subscript  1  refers  to 
conditions  when  corresponding  points  in  A  and  B  are  in  contact  with  each  other. 
Let  Vj(t)  be  the  velocity  of  a  point  A  in  surface  A  that  was  initially  in  contact 
with  B  at  time  tj  and  position  Xj .  Then  the  relative  motion  of  A  may  be  written 
as 


5 


•  .*  *  ,  • « • .  •  .t*  «t*  ./  - 
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■y 

XA  =  X 


J  -  Xj  ♦  /  Vj  dt 


(2.8) 


The  total  distance  travelled  between  tj  and  t  is  R  =  /  |V^|dt 


.  _  \  |vx |  t  >  tj 

10  t  <  t. 


(2.9) 


R  =  0  at  t  =  t, 


As  stated  earlier  the  point  under  consideration  came  into  contact  with  a 
corresponding  point  in  B  at  time  tj.  As  the  point  A  moves  over  the  contact 
surface  a  bond  will  develop.  The  adhesive  strength  will  depend  on  the  extent  of 
time  the  two  points  were  in  contact  with  each  other.  The  time  of  contact,  in 
turn  depends  upon  the  relative  speed  of  the  two  points  at  the  time  of  contact. 

The  adhesive  strength  will  be  proportional  to  the  area  of  the  bond.  Let  e 
represent  the  dimension  of  the  bond  in  the  direction  of  motion.  Before  the  two 
points  come  in  contact  e=0.  After  they  are  in  contact,  they  will  deform  together 
until  the  contact  is  broken.  The  growth  of  this  bond  by  diffusion  may  be 
represented  by  an  evolution  law  as  follows 


e(xj  ,t) 


t  <  t, 


(2.10) 


4s-  <ecTe)  +  — ^ -  (ek~e)  1  l  4 

M  *k 


where  ag,  eQ,  are  constants  that  can  be  chosen  to  represent  a  given 
situation  where  either  the  static  (subscript  s)  or  kinematic  (subscript  k) 
condition  is  expected  to  be  predominant.  The  above  form  is  chosen  to  allow  for 
the  realistic  growth  of  bond  both  at  low  and  high  relative  velocities. 

This  is  illustrated  by  considering  the  steady  state  value  for  e  which  are 
given  by 
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e  «  0 


(2.11) 


For  low  speeds,  (i.e.,  Vj  very  small) 


as  |VJ 

— 8-  »  — —  and  e 

|Vj |  *k 


For  high  speeds 


(2.12) 


The  constants  a{.  and  can  be  used  to  account  for  the  effect  of  tempera¬ 
ture  on  the  rate  of  bond  growth,  i.e.,  if  Tg  is  the  surface  temperature  at  Aj, 
and  Ta  and  are  activation  energies  then  1 


-Vrs, 


ae  =  ae 
s  o 


T./Tc 

*k  "  V  Sl 


(2.13) 


will  account  for  the  rate  of  bond  growth. 

As  relative  motion  begins  and  the  two  points  move  apart,  the  shear  force  at 
the  bond  will  change.  At  first,  the  shear  force  will  increase,  reach  a  maximum, 
and  as  the  bond  continues  to  break,  the  force  will  decrease.  The  above  charac¬ 
teristics  can  be  represented  by  a  function  $  as  follows 
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$  *  $(R/e) 


1 

r(P+i) 


(R/e) 


p  -(R/e) 
e 


(2.14) 


It  may  be  noted  that  the  ratio  R/e  represents  the  relative  length  of  motion  in 
terms  of  bond  length.  This  function  is  shown  in  Fig.  13  for  several  values  of 
the  parameter  p.  R  is  along  the  path  of  motion  and  is  0  by  definition.  Hence 
4  is  finite  if  p  >  -1  which  requires  p  >  -1 . 

Perpendicular  to  the  path,  a  second  weighting  function  is  allowed  to 
represent  the  bond  area.  This  has  been  taken  to  be  (see  Fig.  14) 


W 


e-(n2/2p2) 


(2.15) 


where  rj  is  the  distance  perpendicular  to  the  path.  The  shear  force  acting 
between  the  two  points  is  now  made  proportional  to  the  product  of  $  and  W. 
Further,  the  strength  of  the  bond  must  also  depend  upon  the  normal  stress 
on(xj)  at  Xj  .  Therefore,  the  proposed  nonlocal  law  for  the  shear  force  at  Xj 
is  essentially  complete  with  the  following  representation 


T(X)  -  lis 


/ 


/  4>(R/e)w(n/p)o  (x,) 

-  Eop 


(2.16) 


Recognition  of  the  fact  that  the  frictional  shear  at  X  will  be  in  a  direction 
opposite  to  the  relative  displacement  u  completes  the  required  representation  as 
follows 


t(X)  =  -ys  / 


♦ 

u 


/  ♦(R/e)W(n/p)on(X1)  — 

—  |u| 


dXjdyj 


eoP 


(2.17) 


It  will  be  shown  that  the  constants  present  in  this  formulation  can  be  related 
to  the  static  and  kinetic  coefficients  of  friction.  Consider  the  surface  A 
moving  at  a  constant  velocity  relative  to  B  under  a  constant  normal  stress. 
Then,  it  can  be  shown  that 
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♦  -► 
t(X) 


o 


n 


where  $  is  Che  relative  velocity.  Clearly  the  above  is  a  degeneration  of 
general  representation  to  a  local  law. 


For  slow  speeds  e  ■  e0  and  then 


♦  V 

'00  . 

|v| 

ps  can  be  identified  with  the  static  coefficient  of  friction.  For  high 
velocities  e  «  ek 


and 


t(X)  -  -ii  g 


0 


n 


Hence  the  quantity  u  r*  can  be  identified  with  the  kinetic  coefficient  o 

s  eo 

friction 
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3.  DYNAMICS  OF  STRUCTURAL  SYSTEMS  WITH  INTERRUPTED  INTERFACES 


A  general  structural  system  with  an  arbitrary  number  of  interfaces  is 
considered.  The  friction  forces  developing  at  these  interfaces  during  vibration 
are  treated  in  such  a  way  that  the  influence  of  different  laws  of  friction  may  be 
examined.  The  solution  procedure  leads  to  a  calculation  of  time  history  of  the 
vibratory  motion  either  under  arbitrary  transient  loading  or  as  the  driving 
frequency  is  swept  over  a  given  range  so  that  the  dynamic  response  as  well  as  the 
nature  of  its  stability  can  be  evaluated. 


A.  Formulation 

The  governing  equations  of  motion  for  the  structural  system  may  be 
represented  in  matrix  form  as 


[m]{x}  +  [c]{x}  +  {k]{x}  *  |f}+{F} 


(3.1) 


{ x }  represents  physical  displacements 
{ f }  represents  forces  of  friction 
{f}  represents  external  forces 

[m] ,  [c]  and  [k]  are  the  mass,  damping,  and  stiffness  matrices. 

It  may  be  noted  that  some  of  the  elements  of  { f }  and  { f}  may  be  zero.  Also, 
if  { f }  represent  forces  at  internal  interfaces  only,  the  arithmetic  sum  of  two 
elements  of  { f}  ■  0. 

Let  [*]  represent  the  matrix  of  eigenvectors  {<p}  obtained  from  the  solution 
of  equation  (3.1)  in  which  [c]=0={ f }“{f}  and  [  u>2]  represents  the  corresponding 
eigenvalues . 

Successive  transformations,  as  shown  below,  can  lead  to 


{z}  ♦  [  u.2J{z}  *  {gHG}-[H] 


(3.2) 


in 
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where 


il'V’ 

E*’v' 

hy,v 


W  -  [*lT{y} 

{y}  -  [✓«!{*} 

{g}  -  [*]T[  l//m  ]{f}  (3.3 

{G}  «  [«]T[  l//£  ]  {F} 

[H]  -  [♦]T[c] [♦]{*} 

The  elements  of  {g}  are  functions  of  time  and  the  elements  of  { g}  represent  the 
nature  of  friction  forces. 

The  ith  row  of  Eq.  (2)  can  be  written  as 


zj  +  ujzj  =  g£  ♦  Gj  -  h£  (3.4 

gi  *  g^(z£,z^)  in  view  of  Eq.  (3.3)  represents  the  influence  of  friction  forces 
in  mode  i ; 


G^  ■  G^(t)  represents  the  generalized  force  in  the  ith  mode  and 


h^  =  h^(z^)  represents  the  influence  of  viscous  type  damping  forces  in  the  ith 
mode. 


B.  Solution 

The  solution  to  Eq.  (3.4)  can  be  written  in  terms  of  the  complementary 
solution  and  the  particular  integral,  i.e.,  z  ■  zC  ♦  zP 


ZiC  «  CpVUit  ♦  ci2>e"1Uit 


where  the  constants  C;^  and  C;2^  are  determined  from  the  initial  conditions. 
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The  particular  solution  is  assumed  to  be 


z^  «  Dil>(t>e“*t  ♦  Dp>(t)e'lWi 


-lU^t 


«  D£1)(t)e1<°lt  +  D(i2)(t)e  1U>lt  +  iu)|D(i1)(t)ei“it_Dp)(t)e_i“it| 


As  only  one  function  can  be  arbitrary  in  the  solution  of  Eq.  (3.4),  a 
relationship  between  ^  and  can  be  obtained  by  setting 


D<1>d1“it  ♦  i<2>  _1“:t 


D>-'d  ~  +  D>_'e 


Thus , 


z\  *  id)  |o^1^(t)eiUit  +  D^(t)e 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


•  P 

Differentiating  Eq.  (9)  once  more  to  obtain  and  substituting  in  Eq.  (3.4) 


leads  to  a  solution  for  '  and  . 


z 


iu  •  t 
e  + 


-lu^t 


iu • (t-T ) 
e  1  dr 


(3.10) 


+ 


/  (gi+Gi-Hj 


-iUi(t-T) 
e  1  dT 
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where  gi  -  g^z^r),  z^x)} 

Hi  - 

Gi-G^T) 


It  may  be  noted  that  and  C*^  are  calculated  from  prescribed  initial 

conditions,  i.e.,  {x}  and  { x}  at  t  *  0.  With  {x}  and  {x}  prescribed  at  t  ■  0,  i 
can  be  calculated  from  Eq.  (3);  i.e., 

{ Z1 t=0  "  [*]TC^m] {x}t,0 


{z}t«0  *  { x} t,0 


which  transforms  the  physical  coordinates  to  modal  coordinates. 


C.  Numerical  Application 

The  implicit  integration  procedure  (Eq.  (3.10))  described  above  has  been 
coded  on  a  P.C.  (HP9816)  and  several  example  problems  were  run  to  verify  the 
procedure  and  the  program.  Only  two  examples  are  shown  below  to  illustrate 
certain  aspects  of  the  procedure  including  important  computational  steps  that 
need  special  attention  in  obtaining  reliable  results. 

1 .  Single-Degree-of-Freedom  Example 

The  numerical  solution  of  a  simple  single  degree  of  freedom  system  with  a 
Coulomb  type  of  friction  interface  (Fig.  4)  is  the  first  example.  With  the 
initial  conditions  prescribed  as  x*xQ ,  x»0  at  t*0 ,  the  solution  of  the  governing 
equation  mx  +  kx  *  -f  sgn(x) ,  subjected  to  the  above  initial  conditions  can  be 
shown  to  be 

fx  -  (2N+1 )  — )  cos  u  t  *  (-1)N—  N  <  N. 
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where  N  *  number  of  half  cycles  completed,  truncated  to  the  greatest  integer 
less  than 

o»0t/ir 


Figure  15  shows  the  time  history  of  motion  obtained  by  numerical  integration 

which  is  identical  to  the  exact  solution  at  all  times  until  the  mass  comes  to 
irN 

rest  at  t  *  — —  ~  11  units  of  time  as  shown  in  Fig.  15.  Beyond  this  time,  a  very 
o 

slight  tendency  to  oscillate  is  indicated. 

The  results  from  this  example  confirmed  that  the  code  is  able  to  calculate 
highly  nonlinear  conditions  with  acceptable  accuracy. 

2 .  Two-Degree-of-Freedom  Example 

The  second  example  chosen  is  a  two-degree-of- freedom  system  which  tests  some 
important  aspects  of  frictionally  damped  systems  as  discussed  below. 

Figure  16  is  the  analytical  model  of  a  2D.0.F.  which  has  a  rigid  body  mode 
and  a  flexible  system  mode.  When  mass  2  is  stuck,  the  system  degenerates  into 
the  S.D.O.F.  model  discussed  above.  With  P  *  P0  sin  wt  and  an  assumed  friction 
force  of  pN  •  3P0,  the  system  would  respond  as  a  S.D.O.F.  as  long  as  the 
driving  frequency  is  not  close  to  the  natural  frequency,  /K/m.  At  a  frequency 
close  to  the  natural  frequency,  the  spring  force  would  gradually  build  up  to 


overcome  3PQ  and  mass  2  will  slip.  However,  at  w  *  1.14  /K/m  and  uN  ■  3P0, 
calculations  were  made  using  the  displacement  at  each  time  as  the  criteria  for 
convergence.  The  resulting  time  history  for  the  displacement  of  mass  2  and  the 
friction  force  at  the  interface  are  shown  in  Figs.  17  and  18  respectively.  These 
results  are  unacceptable:  the  mass  2  cannot  begin  to  move  because  the  friction 
force  at  the  interface  has  not  reached  3P0  as  can  be  seen  from  Fig.  18.  This 
is  an  illustration  of  the  influence  of  numerical  analysis  on  the  solution  and 
suggests  the  need  to  exercise  considerable  care  in  the  choice  of  convergence 
criterion. 


\  \  \  *. 
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The  convergence  criterion  was  changed  to  a  check  on  velocity  at  each  step 
and  the  results  show  some  improvements  as  can  be  observed  from  Figs.  19  and  20. 
Mass  2  moves  but  there  is  no  steady  drift.  However  the  frictional  force  is  quite 
erratic.  This  is  an  important  problem  that  points  to  erroneous  results  that 
demand  a  careful  look  into  modeling  friction  dependent  systems. 


Simplified  Nonlocal  Law 


A  final  attempt  was  made  by  using  a  simplfied  version  of  a  nonlocal  law  as 
shown  below 


I F |  <  pN 


IF |  >  pN 


(3.12) 


Convergence  criterion  based  on  velocities  was  used  and  the  results  obtained  are 
shown  in  Figs.  21  and  22.  Clearly,  the  displacement  of  mass  2  is  essentially 
zero  until  the  required  friction  force  exceeds  pN.  Although  the  friction  force 
shows  some  instability,  the  representation  is  superior  to  that  seen  earlier  in 
Fig.  20. 

These  examples  illustrate  the  complexity  of  modeling  friction  at  interfaces, 
even  for  relatively  simple  systems.  Representation  of  the  friction  law  as  well 
as  numerical  methods  have  to  be  chosen  with  utmost  care. 
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I 

i 

'  4.  EXPERIMENTAL  PROGRAM 


|  A.  Configuration 

i  The  photographs  of  the  test  assembly  are  shown  in  Figs.  23  and  24.  Figure 

24  shows  the  details  of  the  test  pieces.  The  entire  friction  assembly  is  mounted 
on  a  slip  plate  connected  to  the  shaker  so  that  predetemined  acceleration  inputs 
may  be  applied  to  the  system.  The  springs  (threaded  rods  connected  to  flexures; 
the  latter  are  cantilever  beams  to  the  tips  of  which  the  threaded  rods  are  con¬ 
nected)  offer  resistance  to  motion  of  the  loaded  mass  and  friction  forces  mani¬ 
fest  at  the  interfaces.  The  normal  load  on  the  block  may  be  varied  by  changing 
the  weight  hanging  from  the  block. 

The  system  had  been  designed  such  that  the  natural  frequency  of  vibration 
for  the  translational  (axial)  mode  of  interest  would  lie  between  100  and  150 
Hz. 


Details  of  the  titanium  test  pieces  are  shown  in  Fig.  25.  The  actual  fric¬ 
tion  surfaces  are  6  mil  hardface  coatings  of  tungsten  carbide,  as  shown. 


B.  Preliminary  Modal  Analysis 

A  preliminary  modal  analysis  was  performed  to  examine  the  dynamics  of  the 
system  after  it  was  assembled  in  place.  The  predominant  mode  of  the  system  is, 
as  designed,  one  of  translation  (see  Fig.  26)  occurring  around  143  Hz,  although  a 
negligible  amount  of  pitching  and  yawing  modes  can  be  present,  especially  when 
the  friction  interfaces  are  disengaged.  There  is  a  "rocking"  mode  of  the  mass  at 
around  73  Hz  with  axial  amplitudes  considerably  less  than  that  at  the  mode  of 
interest.  Damping  in  the  axial  mode  was  predicted  to  be  ~0.003  based  on  a  multi- 
degree-of-freedom  fit  of  the  frequency  response  function. 

A  lateral  yawing  mode  was  found  at  129  Hz  and  could  have  compromised  results 
if  there  had  been  any  significant  asymmetry  in  the  friction  forces.  This  did  not 
prove  to  be  the  case  in  the  subsequent  testing. 

Upon  completing  the  assembly  with  friction  interfaces  in  place,  but  with  an 
external  normal  load  of  26#  on  the  loading  cable,  the  modal  analysis  was  repeated 
and  the  results  are  shown  in  Fig.  27.  The  "rocking"  and  "yawing"  modes  have  been 
suppressed;  the  mode  of  interest  was  found  to  be  almost  critically  damped  in  the 
vicinity  of  148  Hz. 
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C.  Test  Procedure 

Exploratory  sweeps  between  130  and  160  Hz  at  low  acceleration  (<  5g)  were 
made  to  examine  system  response  characteristics.  The  4g  level  was  found  to  be  a 
critical  one  at  all  frequencies  since  the  response  acceleration  waveform  became 
distorted  due  to  the  onset  of  slipping  conditions.  The  input  level  was  found  to 
be  the  more  sensitive  parameter  and  hence  a  procedure  was  selected  which  required 
the  input  level  to  be  increased  gradually  in  steps  while  the  frequency  was  held 
constant.  This  was  performed  at  four  frequencies  namely,  135,  143,  150  and  155 
Hz.  The  input  level  was  varied  from  the  level  at  which  slipping  was  initiated, 
until  either  a  divergent  condition  or  10  g  occurred.  The  input  acceleration, 
moving-block  absolute  acceleration  and  flexure  strain  gage  output  were  recorded 
on  a  visicorder.  A  separate  static  test  was  performed  to  determine  the  constant 
of  proportionality  to  convert  the  strain  gage  conditioner  output  to  relative 
displacement  of  the  moving  block.  The  constant  was  found  to  be  approximately  1 
volt  per  0.001  inch. 


D.  Test  Results 

Typical  traces  of  the  input  and  responses  for  slipping  conditions  are  given 
in  Figs.  28  and  29.  Figure  28  shows  the  responses  at  135  Hz/6g  input  and  Fig.  29 
shows  the  response  for  143  Hz  input  frequency  at  5g  and  8g  levels.  At  the  8g 
level  a  stable  and  an  unstable  condition  (which  occurred  without  warning)  are 
shown. 

The  friction  surfaces  were  characterized  before  and  after  the  series  of 
tests  by  means  of  a  Sloan  Dektak  II  prof ilometer .  Figure  30  shows  the  asperity 
profiles  before  the  test  series  in  the  x  and  y  directions,  and  Fig.  31  the 
changed  profiles  taken  after  the  tests.  The  before-test  profiles  show  an 
important  component  of  80  micron  scale  with  an  amplitude  of  13  to  18  microns, 
apparently  isotropic  in  the  x,  y  coordinates.  After  test,  the  80  micron  scale 
asperities  are  still  visible  but  of  much  smaller  amplitude,  and  with  a  smaller 
amplitude  in  the  x  direction  (the  direction  of  sliding). 

The  response  at  resonant  frequency  to  an  input  displacement  of  Xo  =  .0002  ft 
in  amplitude  (5g  equivalent  acceleration)  is  shown  in  Figures  2,  3,  32  and  33 
where  these  responses  are  shown  in  dimensionless  form.  Since  the  term  Mw2/K  in 
equation  2.3  becomes  equal  to  1  at  a  resonance  with  small  damping,  x  ■  q  Xo. 

From  Figure  2,  q  ■  .15;  therefore  x  becomes  .00003  ft  or  .00036  inches,  comparing 
well  with  the  measured  value  of  .00023  inches  in  Table  1  (5g/143  Hz).  This  is  a 
good  agreement  considering  the  critical  nature  of  the  response  at  the  undamped 
resonant  frequency  of  the  system. 
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At  resonance j^Che  value  of  xM/KXo  *  1  as  can  be  seen  from  Figure  3.  This 
can  be  solved  for  x  leading  Co  a  value  of  5g  for  Che  absoluce  acceleration, 
comparing  closely  wich  Che  6g  in  Table  1.  The  response  was  also  compuCed  for 
equivalenC  4g  and  8g  input  accelerations.  In  Che  4g  case  no  slipping  occurred 
and  in  Che  8g  case  divergence  of  Che  response  occurred.  This  compares  favorably 
with  Che  results  obtained  from  tests  (see  Table  1). 
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5.  CONCLUSIONS 


The  following  observations  summarize  the  principal  findings  from  the  tests 
as  well  as  from  a  comparison  between  test  data  and  corresponding  analysis.  These 
observations  are  to  be  understood  as  valid  for  the  frequency  and  g  level  ranges 
tested.  It  must  be  noted  that  the  measured  acceleration  is  absolute  and  is 
compared  with  corresponding  calculation.  Relative  acceleration  from  the  calcula¬ 
tions  is  also  shown. 

1)  Slipping  at  the  interfaces  occurs  at  and  above  a  certain  g  level 
independent  of  the  excitation  frequency  (see  Table  1). 

2)  Above  a  certain  g  level  the  amplitudes  of  relative  slip  motion  experi¬ 
enced  divergence  limited  entirely  by  nonlinearities  due  to  the  increase 
in  flexure  stiffness  at  large  amplitude. 

3)  As  the  mass  undergoes  sticking  and  slipping  the  resulting  accelerations 
are  sudden  and  large.  This  creates  an  impulse  type  of  input  which  in 
turn  induces  higher  frequency  responses  seen  as  ringing  (See  Fig.  28  and 
29). 

4)  Responses  can  be  nearly  sinusoidal  above  certain  g  levels. 

5)  Rubbing  was  observed  to  have  occurred  only  on  a  portion  of  the  total  area 
of  contact.  It  is  likely  that  the  engaging  surfaces  were  not  perfectly 
aligned  due  to  machining  tolerance. 

6)  As  the  input  displacement  reaches  a  maximum  (Fig.  28,  32),  provided  the  g 
level  is  sufficiently  high,  slipping  begins  to  occur  and  reaches  a 
maximum  somewhat  later  than  the  time  at  which  the  input  acceleration 
reaches  its  maximum. 

7)  As  the  input  displacement  is  returning  to  zero  upon  completing  the  first 
half  cycle,  the  relative  motion  at  the  interfaces  remains  constant 
(i.e.  the  masses  are  stuck). 

8)  As  the  input  displacement  continues  to  increase  beyond  the  first  half 
cycle,  the  masses  begin  to  experience  relative  slip  again  and  try  to 
follow  the  input,  although  lagging. 

9)  At  the  instant  of  sticking,  the  velocity  suddenly  drops  to  zero  (Fig.  2) 
from  a  finite  value,  resulting  in  a  sudden  increase  in  acceleration. 

This  is  shown  by  large  discontinuities  in  the  acceleration  traces.  Where 
the  discontinuities  appear  several  times,  the  interpretation  is  that 
local  stick-slips  continue  until  a  full  slip  occurs. 
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■Ill 

■III 

■III 

■III 

■■■■assail 

■■■■ 

■■■■ 

■HI 

■HI 

■■■■ 

Overall  Arrangement 


CONNECTING  ROD 


Friction  Test  Assembly  Showing  Test  Pieces 
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DRILL  THROUGH  2  HOLES 
#8  (0  199  DIA) 
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a)  UPPER  TEST  PIECE 


DRILL  THROUGH  4  HOLES  #8  (0  199D) 
82 ’  COUNTERSINK  TO  0  42  MAX  DIA 


>  THESE  SURFACES  COATED  WITH 
0  006  THICK  TUNGSTEN  CARBIDE 
HARDFACE 
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DRILL  #38  (0  1015) 
0  18  DEEP 
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b)  LOWER  TEST  PIECE 


Fig.  25  Details  of  Titanium  8-1-1  Test  Pieces 
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